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The exact joint muhifractal distribution for the scaling and winding of the electrostatic potential 
lines near any conformally invariant scaling curve is derived in two dimensions. Its spectrum /(a, A) 
gives the Hausdorff dimension of the points where the potential scales with distance r as i/ ~ r" 
while the curve logarithmically spirals with a rotation angle = Alnr. It obeys the scaling law 
/(a, A) = (1 + \^)f (q) - with Q = a/(l + A^) and 6 = (25 - c)/12, and where f{a) = f{a, 0) 
is the pure harmonic measure spectrum, and c the conformal central charge. The results apply to 
0{N) and Potts models, as well as to SLE^. 
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The geometric description of the random fractals aris- 
ing in Nature is a fascinating subject. Among these, 
the study of the particular class of random clusters or 
fractal curves arising in critical phenomena has led to 
fundamental advances in mathematical physics. In two 
dimensions (2D), conformal field theory (CFT) has in 
particular demonstrated that statistical systems at their 
critical point produce conformally invariant (CI) fractal 
structures, examples of vi^hich are the continuum scaling 
limits of random walks (RW), i.e., Brownian motion, self- 
avoiding walks (SAW), and critical Ising or Potts clusters. 
A wealth of exact methods has been devised for their 
study: Coulomb gas, conformal invariance, and quantum 
gravity methods |, |, |, |] . Recently, rigorous meth- 
ods have also been developed, with the introduction of 
the Stochastic Lowner Evolution (SLE) process, which 
mimics the wandering of critical cluster boundaries, and 
gives a probabilistic means of studying their random ge- 
ometry (H. 

A refined way of accessing this random geometry is pro- 
vided by classical potential theory of electrostatic or dif- 
fusion field near such random fractal boundaries, whose 
self-similarity is reflected in a multifractal (MF) spec- 
trum describing the singularities of the potential, also 
called the harmonic measure. In 2D, the first exact ex- 
amples appeared for the universality class of random or 
self- avoiding walks, and percolation clusters, which all 
possess the same harmonic MF spectrum [p| (see also [Q), 
in contradistinction to higher dimensions |^| . The general 
solution for the potential distribution near any CI frac- 
tal in 2D, obtained in [^, depends only on the so-called 
central charge c, the parameter labeling the universal- 
ity class of the underlying CFT (see also [ll). This 
solution can be generalized to higher multifractal corre- 
lations, like the joint distribution of potential on both 
sides of a simple scaling path . 

The important question remains of the geometry of 
the equipotential lines near a random (CI) fractal curve. 



They are expected to wildly rotate, or wind, in a spi- 
ralling motion which closely follows the boundary it- 
self. The key geometrical object is the logarithmic spi- 
ral, which is conformally invariant. The MF description 
should generalize to a mixed multifractal spectrum, ac- 
counting for both scaling and winding of the equipoten- 
tials 

In this Letter, we obtain the exact solution to this 
mixed MF spectrum for any random CI curve. In partic- 
ular, it is shown to be related by a scaling law to the usual 
harmonic MF spectrum. We use conformal tools (fusing 
quantum gravity and Coulomb gas methods), which al- 
low the description of Brownian paths interacting and 
winding with CI curves, thereby providing a probabilis- 
tic description of the potential map. 

Harmonic Measure and Rotations. Consider a single 
(CI) critical random cluster, generically called C. Let 
H (z) be the potential at the exterior point z G C, with 
Dirichlet boundary conditions H {w E dC) = on the 
outer (simply connected) boundary dC of C, and H{w) = 
1 on a circle "at oo" , i.e., of a large radius scaling like the 
average size R of C. As is well-known, H (z) is identical 
to the probability that a Brownian path starting at z 
escapes to "oo" without having hit C. 

Let us now consider the degree with which the curves 
wind in the complex plane about point w and call (p{z) = 
arg (z — w). The multifractal formalism 0, |l5|, [l^ , 0, 
here generalized to take into account rotations |l3|, char- 
acterizes subsets dCa,\ of boundary sites by a Holder ex- 
ponent a, and a rotation rate A, such that their potential 
lines respectively scale and logarithmic spiral as 

(z -> u; e aCa,A) ~ r", 
(p{z w E dCa,\) ~ Aln r , 

in the scaling limit oq <^ r — \z — w\ R, where ag 
is the lattice mesh, if any. The Hausdorff dimension 
dim {dCa,\) = f (a, A) defines the mixed MF spectrum. 
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which is CI since under a conformal map both a and A 
are locally invariant. 

Reversing the escaping Brownian path which repre- 
sents the potential, one can also consider the harmonic 
measure H{w,r), which is the probability that such a 
path starting at distance R first hits the boundary in the 
disk B{w, r) of radius r centered at w; € dC, and fi{w, r) 
the associated winding angle of the path down to distance 
r from w. The mixed moments of H and e*^, averaged 
over all realizations of C, are defined as 



Zn,p=l J2 ^"KOexp(p^(w;,r))\ «(r/i?) 

(1) 

where the sum runs over the centers of a covering of the 
boundary by disks of radius r, and where n and p are real 
numbers. The scaling limit involves multifractal scaling 
exponents t (n, p) which vary in a non- linear way with 
n and p [|l3[ Q |l5[ [l^, They obey the symmetric 
double Legendre transform 

a = ^{n,p), X-^in,p), 
/(a, A) — an + Xp — T {n,p) , 

n = f^(-,A), p=|{(a,A). (2) 

Because of the ensemble average values of f{a,X) 
can become negative for some domains of a, A. 

Exact Mixed Multifractal Spectra. Each 2D confor- 
mally invariant random statistical system can be labelled 
by its central charge c, c < 1. Our main result is the fol- 
lowing exact scaling law: 



/(a, A) 
b 
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1 + A2 



bX' 



(3) 
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where / (a) = / (a, A = 0) is the usual harmonic MF 
spectrum in the absence of prescribed winding, first ob- 
tained in PI, which can be recast as: 



/(«) 



ba^ 



2a - 1 

We thus arrive at the very simple formula: 
f{a,X)=a+' 



2a - 1 - A2 



(4) 



(5) 



Notice that by conformal symmetry sup;^/(a,A) = 
/(a, A = 0), i.e., the most likely situation in the absence 
of prescribed rotation is the same as A = 0, i.e. winding- 
free. The domain of definition of the usual /(a) (^) is 
a > 1/2 [|, 18 1, thus for A-spiraUing points Eq. (g) gives 




FIG. 1: Universal multifractal spectrum f{a, X) for c — 
(Brownian frontier, percolation EP and SAW), and for three 
different values of the spiralling rate A. 



in agreement with a theorem by Beurling ||l^, [l^l . 

There is a geometrical meaning to the exponent a. For 
an angle with opening 9, a = n/O, thus the quantity 7r/a 
can be regarded as a local generalized angle with respect 
to the harmonic measure. The geometrical MF spectrum 
of the boundary subset with such opening angle 9 and 
spiralling rate A reads from (||) 

f{9,X)^f{a = -,X) = ^ + b-b-i- 



2tt 



As in (^), the domain of definition in the 9 variable is 
< 9 < 6'(A), with 61(A) = 27r/(l + A^). The maximum 
is reached when the two frontier strands about point w 
locally collapse into a single A-spiral, whose inner opening 
angle is 9{X) 

In the absence of prescribed winding (A = 0), the max- 
imum Dep = i?Ep(0) = sup„/(a,A — 0) gives the di- 
mension of the external perimeter of the fractal cluster, 
which is a simple curve without double points, and may 
differ from the full hull P, p^ . Its dimension reads |^ 
Z^EP = ^(l+&) — ^■\/6(fe — 2). This corresponds to typical 
values d — a{n = 0,p = 0) and 9 = 7r/d = 7r(3 — 2Z?ep)- 

For spirals, the maximum value Dep{X) = snp^f{a, A) 
still corresponds in the Legendre transform (^ to n 
0, and gives the dimension of the subset of the external 
perimeter made of logarithmic spirals of type X. Owing 
to (PI) we immediately get 



^Ep(A) = il + X'')DEP~bX^ 



(7) 



a>-(l + A2) 



(6) 



This corresponds to scaled typical values d(A) = (1 + 
A2)d, and ^(A) = 9/{l + A^). Since 6 > 2 and i:»EP < 
3/2, the EP dimension decreases with spiralling rate, in 
a simple parabolic way. 

Fig. 1 displays typical multifractal functions /(a, A ; c). 
The example choosen, c — 0, corresponds to the cases 
of a SAW, or of a percolation EP, the scaling limits of 
which both coincide with the Brownian frontier 
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FIG. 2: Left-sided multifractal spectra f{a, A) for the limit 
case c — 1 (frontier of a Q = 4 Potts cluster or SLEk=4). 



The original singularity at a = ^ in the rotation free 
MF functions /(a,0), which describes boundary points 
with a needle local geometry, is shifted for A 7^ towards 
the minimal value (^. The right branch of / (a, A) has 
a linear asymptote lim^^+oo / (a, A) /a = —(1 — c)/24. 
Thus the A-curves all become parallel for a +00, i.e., 

— > 0"*", corresponding to deep fjords where winding is 
easiest. 

Limit multifractal spectra are obtained for c — 1, 
which exhibit exact examples of left-sided MF spectra, 
with a horizontal asymptote f (a ^ +00, A ; c = 1) = 

1 — iA^ (Fig- 2). This corresponds to the frontier of 
a Q = 4 Potts cluster (i.e., the SLE„-4), a universal ran- 
dom scaling curve, with the maximum value Z^ep = 3/2, 
and a vanishing typical opening angle = 0, i.e., the "ul- 
timate Norway" where the EP is dominated by "fjords" 
everywhere [|[ p^ . 

Fig. 3 displays the dimension _Dep(A) as a function 
of the rotation rate A, for various values of c < 1, cor- 
responding to different statistical systems. Again, the 
c = 1 case shows the least decay with A, as expected 
from the predominence of fjords there. 

Conformal Invariance and Quantum Gravity. We now 
give the main lines of the derivation of exponents t (n,p), 
hence /(a. A), by generalized conformal invariance. By 
definition of the _ff-measure, n independent Brownian 
paths B, starting a small distance r away from a point 
w of the frontier dC, and diffusing without hitting dC, 
give a geometric representation of the n*'' moment, i?", 
in Eq. (|^) for n integer. Convexity yields analytic continu- 
ation for arbitrary rt's. Let us introduce an abstract (con- 
formal) field operator $acAn characterizing the presence 
of a vertex where n such Brownian paths and the clus- 
ter's frontier diffuse away from each other in a mutually- 
avoiding configuration noted dC An ; to this operator 
is associated a scaling dimension x{n). To measure rota- 
tions as in moments (|^) we have to consider expectation 
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FIG. 3: Dimensions _Dep(A) of the external frontiers as a 
function of rotation rate. The curves are indexed by the 
central charge c, and correspond respectively to: loop-erased 
RW (c = — 2;SLE2); Brownian or percolation external fron- 
tiers, and self-avoiding walk (c — 0;SLEg/3); Ising clusters 
(c = i ; SLE3); (5 = 4 Potts clusters (c = 1 ; SLE4). 



values with insertion of the mixed operator 

<facA„eP--s(^'^^") ^ a- ^ ^ (^^p) + 2, (8) 

where arg(9C A n) is the winding angle common to the 
frontier and to the Brownian paths, and where x{n,p) is 
the scaling dimension. One has x(n,p = 0) = a;(n), and 
T (n,p = 0) = r (n) = x (n) — 2. 

Let us now use a fundamental mapping of the CFT 
in the plane to the CFT on a fluctuating abstract 
random Riemann surface, i.e., in presence of 2D quantum 
gravity (QG) [po| . Two universal functions U and 
acting on scaling dimensions, describe this map: 



U{x) = 



■ 7 



1 



7 



1 

V(x) = - — 
^ ' 4 1 



7 



■7 



(9) 



with F (a;) = ?7 (i (a; + 7)) [|, ||. The parameter 7 is the 
solution of c = 1 — 67^(1 — 7)"^, 7 < 0. 

For the purely harmonic exponents x{n), describing 
the mutually-avoiding set dC A n, we have [p[ ^ 

x{n) = 2V [2U-'^ {£1) + U-^ (n)] , (10) 

where (x) is the positive inverse of U 

2C/-1 (x) = ^4(1 -7)a; + 72 + 7 . 

In ([To|), the arguments xi and n are respectively the 
boundary scaling dimensions (b.s.d.) of the simple path 
Si representing a semi-infinite random frontier (such that 
dC = Si A Si) , and of the packet of n Brownian paths, 
both diffusing into the upper half-plane H. The func- 
tion U^^ maps these half-plane b.s.d's to the correspond- 
ing b.s.d. 's in quantum gravity, the linear combination 
of which gives, still in QG, the b.s.d. of the mutually- 
avoiding set dC An = (AiSi)^ A n. The function V finally 
maps the latter b.s.d. into the scaling dimension in R^. 
The path b.s.d. afi obeys [/-^ (xi) = (1 - 7)/2 M. 
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It is now useful to consider k semi-infinite random 
paths iSi, joined at a single vertex in a mutually-avoiding 

k 



siar configuration 5fe =<Si A iSi A • • • iSi= (ASi)''. Its 
scaling dimension can be obtained from the same b.s.d. 
additivity rule in quantum gravity, as in ( p^ ) ^ 

x{Sk) = 2V[kU-'{x\)] . (II) 

The scaling dimensions ([l0| ) and ( [TT| ) coincide when 



(n) 



(12) 



a;(n) = a;(5fc(„)), fc(n) = 2 + — — — 
^ ' (7 ^ (xi) 

Thus we state the scaling star- equivalence dC A n -4=> 
'5fc(„), o/ two simple paths Si avoiding n Brownian mo- 
tions to k(n) simple paths in a mutually-avoiding star 
configuration, an equivalence which will also play a es- 
sential role in the complete rotation spectrum (^). 

Rotation scaling exponents. The Gaussian distribution 
of the winding angle about the extremity of a scaling 
path, like 5i, was derived in using exact Coulomb 
gas methods. The argument can be generalized to the 
winding angle of a star Sk about its center |2^ , where one 
finds that the angular variance is reduced by a factor 
(see also ^^). The scaling dimension associated with 
the rotation scahng operator ^^^e^^''^^'^'') is found by 
analytic continuation of the Fourier transforms evaluated 
there I2I: 



x{Sk;p) = x(Sk) 



2 



i.e., is given by a quadratic shift in the star scaling ex- 
ponent. To calculate the scaling dimension (|^), it is suf- 
ficient to use the star-equivalence ( |l2|) above to conclude 
that 



x{n,p) = x(S^.n)\p) = x{n) 



P 



I — 7 k^(n) ' 



which is the key to our problem. Using Eqs (|12|), (10) 
and (^) gives the useful identity: 

1(1- j)k^{n) ^x{n) -2 + b , 



where /(a) = an — T{n), with a = dT(n)/dn, is the 
purely harmonic MF function. It depends on the natural 
reduced variable a d la Beurling (a G [5, +00)) 



a 



dx , ^ I I 
A2 = d;^(")=2 + 2 



2n + 6 - 2 ' 



whose expression is found explicitly from (p^). Whence 
Eq.(|), QED. 

0{N) and Potts models, SLE^. Our results apply to 
the critical 0{N) loop model, or to the EP's of critical 
Fortuin-Kasteleyn (FK) clusters in the Q-Potts model, 
all described in terms of Coulomb gas with some coupling 
constant g . SLE^ paths also describe cluster frontiers 
or hulls. One has the correspondence k — 4/g, with 
a central charge c = (3 — 2g)(3 — 2g') = i(6 — k)(6 — 
k'), symmetric under the duality gg' = I or kk' — 16. 
This duality gives FK-EP's as some simple random 0{N) 
loops, or, equivalently, the SLEk'<4 as the simple frontier 
of the SLE,>4 § @. 
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as is a careful reading of the manuscript by T. C. Halsey. 



witlLii 



[8] 

[9] 
[10] 

[11] 
il2] 



see, e.g., B. Nienhuis, in Phase Transitions and Criti- 
cal Phenomena, edited by C. Domb and J. L. Lebowitz, 
(Academic, London, 1987), Vol. 11. 

J. C. Cardy, Lectures on Conformal Invariance and Per- 
colation, Chuo University, (2001). 

B. Duplantier and I.K. Kostov, Nucl. Phys. B340, 491 
(1990). 

B. Duplantier, Phys. Rev. Lett. 81, 5489 (1998). 
O. Schramm, Israel Jour. Math., 118, 221 (2000). 
B. Duplantier, Phys. Rev. Lett. 82, 880, 3940 (1999). 
G.F. Lawler and W. Werner, Ann. Probab. 27, 1601 
(1999) ; J. Eur. Math. See. 2, 291 (2000) ; J. Cardy, 
J. Phys. A 32, L177 (1999). 

M.E. Gates and T.A. Witten, Phys. Rev. A 35, 1809 
(1987). 

B. Duplantier, Phys. Rev. Lett. 84, 1363 (2000). 

G.F. Lawler, O. Schramm and W. Werner, Acta Math. 

187, 275 (2001). 

M. B. Hastings, Phys. Rev. Lett. 88, 055506 (20021 

B. Duplantie r, J. Stat. Phys., to appear, cond 



_ 1 (2-7)" _ 25- 



"2 — T^^T 

multifractal result: 



T3- 



Rpcalling 



we arrive at thp 



mat/0207743 



I 



T{n,p) = T(n) - - 
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where T(n) = x{n) — 2 corresponds to the purely har- 
monic spectrum with no prescribed rotation. 

Legendre transform. The structure of the full r- 
function (^3|) leads by a formal Legendre transform ^ 
directly to the identity 

/(a. A) = (I + A2)/(a)^6A2 , 
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